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Abstract 

A parametrization of multipartite separable states in a finite-dimensional 
Hilbert space is suggested. It is proved to be a diffeomorphism between 
the set of zero-trace operators and the interior of the set of separable 
density operators. The result is applicable to any tensor product de- 
composition of the state space. An analytical criterion for separability 
of density operators is established in terms of the boundedness of a 
sequence of operators. 



Introduction 

We explore multipartite quantum systems with state space eg) • • -<g>H( N \ 
assuming the factors are Hilbert spaces of finite dimension: 

dim# (1) ,...,dim# (Ar) < oo 

A pure state of the system is a state whose density operator is one- dimensional 
projector. Among pure states we there states of special kind: 

P = p( 1 )®...®pW 
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where P^ k > is a one-dimensional projector in H^ k \ A density operator p 
representable as a convex combination 



where P is the set of all pure product states, is said to be separable. Such 
description of the set of separable density operators is not constructive as 
their appropriate representation is not unique. In this paper a uniform rep- 
resentation of separable density operators p in the form 



is suggested. This representations turns out to be unique, therefore, it can be 
treated as coordinatization of the set of separable density operators. Along 
these lines, a criterion to determine the separability of a given density oper- 
ator p is formulated. 

The paper is organized as follows. In Section [1] separable density oper- 
ators are treated as barycenters of continuously distributed unit masses on 
the set P of pure product states. This gives rise to the equation (J5]) on which 
the coordinatization of separable density operators is based. In Section [2] a 
scalar function (jH]) is studied discriminating the separability in many-particle 
setting. This makes it possible to formulate a criterion of multipartite sepa- 
rability in terms of boundedness of certain sequence of operators. In Section 
[3] the range of the suggested coordinatization of separable density operators 
is determined. The range turns out to be the interior of the set of separable 
density operators, and the coordinatization is proved to be a diffeomorphism. 

1 Separable density operators viewed as barycen- 
ters 

In this section we describe how density operators in a finite-dimensional 
Hilbert space can be described as barycenters of continuous probability dis- 
tributions on a set of one- dimensional projectors. 

1.1 From finite sums to continuous distributions 

We study a multipartite quantum system whose state space H is a finite 
tensor product of finite-dimensional Hilbert spaces H = <g>- ■ ■ <g)H( N '. A 




p 
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density operator p in H is called separable whenever it is an element of the 
convex hull of the set P of pure product states. Since P is the set of extreme 
points of the set of all separable density operators, Choquet theorem makes 
it possible to write down any such p as an integral 

P = ^PdMP) (i) 

over certain probability measure /i on P. In other words, we may repre- 
sent p as the bary center of certain distribution of unit mass on the set P. 
Such representation is known to be essentially non-unique. According to 
Caratheodory theorem this mass may have discrete finite distribution 

p = ^2w a p a 

with the measure p = ^2w a Sp a being a sum of atomic measures 6p a . 

We emphasize that representation ([T]) also comprises continuous distri- 
butions 

p = J P w (P)da(P) (2) 

where w(P) is a positive continuous function and the integration is carried 
out over the probability measure da invariant with respect to all local unitary 
transformations in H. 

So, a density operator p in the product space H = g) • • • <g) is 
separable when it is a barycenter of a continuous distribution w(P) on the 
set P of pure product vectors. In other words, that means that a function 
w(P) exists satisfying the equation 

J Pw(P)da(P) - p = (3) 
1.2 Main equation 

In this section we confine ourselves to distributions w(P) of specific form and 
explore the existence of appropriate solutions. 

Let us first try to find solutions of ([3]) of the form 

W (P) = e Tr ( sp ) (4) 
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where B is an Hermitian operator. Then solving the equation reduces to 
finding an operator X satisfying what we call main equation 

J Pe Tr ( xp W(P) - p = (5) 

In order to verify the existence of a solution of (jHJ), introduce the function 

G(X) = J e Tr(xp W(P) - Tr(Xp) (6) 

whose gradient is 

VG{X) = J Pe Tr(xp W(P) - p (7) 

So, solving (jSJ) reduces to finding an extremum of the function G(X). Taking 
into account that G(X) is convex, the solution of (jSJ) exists only when the 
minimum of G(X) exists. In particular, when p is entangled, there is no way 
for the function G(X) to have minimum: 

Proposition 1. If the density operator p is entangled, then 

MG(X) = -oo 

Proof. The set of all separable states is closed, therefore, if p is not separable, 
there exists a hyperplane, defined by a self-adjoint operator X such that 
VP G P Tr(PX) < 0, while Tr(pX) > 0. Denote 

a = Tr(pX) 

b = maxTr(PX) 

Then a > 0, b < 0, so e Tr ( xp ) < e kb and 

G{kX) < J e kb da(P) - ka -> -oo 

as k — > oo. □ 

So far, finding conditions for the existence of minimum becomes essential 
to judge if p is entangled or not. 
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2 Exponential distributions and their approx- 
imations 



2.1 An interlude on minima of convex functions 

Proposition 2. Let a convex function F(X) has minimum on a finite- 
dimensional space £ and the minimal point is unique. Then 

lim F(X) = +oo (8) 

Proof. With no loss of generality assume the minimum to be attained at 0, 
then F(X) > F(0) for any X ^ 0. Suppose lim X ->oo F(X) ^ +oo, then 
there exists such M and such sequence Xk — >■ oo that F(Xk) < M. Then 

ll^fcll 

Denoting 

\\^k\\ 

we obtain a bounded sequence in the finite-dimensional space S, which con- 
tains a converging subsequence. With no loss of generality denote this sub- 
sequence Ej. and its limit 

E = \imE k , ||^|| = 1 

Then 

Ek = AfcX fc = (1 — Afc) + XkX k 
so, due to the convexity of F 

F(E k ) < (l-Xk)F(0) + X k F(Xk) 

therefore 

F(E k ) < (l-\k)F(0) + X k M 

for sufficiently large k. Taking the limit and recalling that — > 0, and using 
the continuity of convex functions (see, e.g. pQ), we come to the contradic- 
tion: F(E) < F(0). □ 
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So far, we obtained a necessary condition for the unique minimum of a 
convex function to exist. This condition (jSJ) is also sufficient for a minimum 
(not necessarily strict) to exit. In the sequel we shall need a more verifiable 
sufficient condition: 

Proposition 3. Let F(X) be a convex function on a finite- dimensional space 
£. If 

MX ± lim F(tX) = +00 (9) 

t— >oo 

then there exists minimum of F(X) on £. 



Proof. First prove that 



lim F(X) = +00 

X— >oo 



Suppose this is not the case, then there exists a number M and a sequence 
Xk — > 00 such that F(Xk) < M for all k. Denoting 

„ Xu 
z? 



\\Xk\\ 

we obtain a bounded sequence in the finite-dimensional space S, which con- 
tains a converging subsequence. With no loss of generality denote this sub- 
sequence Ek and its limit 

E = \imE k , ||^|| = 1 

Now take an arbitrary t > and N such that ||Xfc|| > t for all k > N. Then 
H^fcll = t < \\Xk\\. Denote 

then < Afc < 1 and tEk = XkXk- We may treat tEk as convex combination 

tEk = (1 — Afc) + AfcXfc 

and apply Jensen's inequality 

F{tE k ) < (1 - A fc )F(0) + A fc F(X fc ) < F(0) + M 

for all k > N. Taking the limit, we get F(tE) < F(0)+M, which contradicts 
with the statement F(tE) — > 00. 
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Now let us prove the existence of minimum. For that, consider the set 



K={X\ F(X) < F(0)} 

The function F is convex (and hence continuous), so the set K is closed. As 
proved above, limx-s-oo F(X) = +00, therefore the set K is bounded. 

The space £ is finite-dimensional, the set K is compact, therefore the 
function F always attains minimum on K, which is its minimum on the 
whole space £ . □ 

So far, we have obtained a sufficient condition for the existence of the 
minimum of a convex function. 

2.2 Multipartite setting 

Now return to our initial setting. We are dealing with integrals over the set P 
of pure product projectors in a tensor product space H = ® • • • ® H^ N \ 
We shall need the following: 

Proposition 4 (Multipolarization identity). Any quadratic form Y in 
H is completely defined by its values on P. 

(e* 1 * ® ...® e^lKl/* 1 ) ® • • • ® /W) = 

3 IN N 

fc(i),...fcCJV)=o \/=o /=o 

(10) 

Proof. Verified by direct calculation. □ 

Remark. In the sequel we shall use the following consequence of this propo- 
sition: 

VP G P Tr(XP) = X = (11) 

Now return to the convex function (El) 



G(X) = J e Tr(xp) da(P) - Tr(Xp) 



Proposition 5. When the function G(X) has minimum, this minimum is 
strict. 
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Proof. Suppose there are two minimal points B ^ B\ of G(X). Denote 
V = Bi — B and consider a family B t = B + Vt. Consider the function 

g(t) = G(B t ) = J e Tr ( BtP W(P) - Tr(i^) 
p 

for t G [0, 1]. This function is constant because G is convex. Therefore its 
second derivative vanishes, but: 

g"(t) = J e Tr(BtP) (Tr(V r P)) 2 d ( x(P) (12) 
p 

In the meantime V ^ 0, therefore J p {Tr(VP)) 2 da{P) > due to (JTTJ), so 

g"{t) > — contradiction. □ 

However, the existence of the minimum of G(X) still can not be directly 
verified. In the meantime, the approximations of G(X) by the family of 
convex functions 

G k (X) = jf (l + : ^^) 2 'W) - Tr(Xp) (13) 
possess the following property: 

Proposition 6. Whatever be p, for any k there exists minimum of the func- 
tion Gt{X) and this minimum is strict. 

Proof. First prove that the minimum exists. Fix an arbitrary I / and 
apply the sufficient condition (j9]) by showing that GkitX) — > +oo. Let 

2k 

then 

G k {tX) = J Ti(Y P) 2k da(P) - t Tr(XA) 
For all F^O, introduce 



Jp Tr(YP) 2k da(P) 
Tr(F 2 



N(Y) = JP > 
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Being homogeneous, N(Y) is completely denned by its values on the compact 
set defined by Tr(F 2 ) = 1. Let us prove that N(Y) is strictly positive. 
Suppose N(Y) = for some Y, then Tt(YP) = for all product one- 
dimensional projectors P G P. That means, for all vectors . . . , 

(e (1) ® • ■ ■ <g> e {N )\Y\e m g> • • • ® e (Af) ) = 

Then, by virtue of (fTTj), F = 0. 

Being continuous function defined on a compact set, N(Y) attains its 
minimal value, denote it a. Then 

k 



G k {tX) > aTr(Y 2 ) k -t Tr(Xp) = at 2h (tt (t~ 1 I+ j -t Tr(Xp) 

since a > as N(Y) was proved to be strictly positive. So the minimum of 
Gk(X) exists. 

In order prove that the minimum is strict, we proceed in a way similar to 
Proposition O with the only difference that the function g(t) has the form: 

git) = G k (B t ) = J ( 1 + : ^|^) 2 "MP) " Tr(^) 
p 

and we check its 2k-s derivative (rather than the second one) 

9 m (t) = J |^ (Tr(VP)) 2 *d<x(P) 
p 

and obtain the same contradiction. □ 
Corollary. Whatever (separable or entangled) be p, it decomposes into 



-> +oo 



J Pw(P)da(P) (14) 



Trf-BPr _ 



with 

»( p ) = ( 1+ lrj (15) 

where B is the minimal point of the function G^, that is, VGk{B) = 0. 

We emphasize that the obtained decomposition f fT4|) is in general not 
barycentric as the density ffl5l) may take negative values. Furthermore, for 
entangled p the density ( fTBT) will always take both positive and negative 
values. Let us consider it in more detail. 



9 



2.3 The convergence of the approximations 

We start from the main equation (jSJ) 

Pe Tr ( xp W(P) 



P 







which may not have solutions and replace it by a sequence of its binomial 
approximations 

Tr(XP)\ 2/c_ 



P I 



2k 



dcr(P) 



P 







(16) 



each always having a unique solution. So, we have to study the conditions 
when these approximations turn to the main equation. For that, we ex- 
plore the convergence Gk — > G of functions associated with the equations in 
question on the set £ of self-adjoint operators in H. 

Theorem 7. If the function G(X) has strict minimum on £ attained in B, 
then the sequence Bk of the minimal points of Gk{X) converge to B. 

Proof. First prove that the convergence Gk — > G is uniform on any compact 
subset of £. For any < a < n direct calculation yields 



max 

\x\<a 

therefore 



f \ n 
e* - ( 1 + - 

n 



max 



1 + - 



i-r } 



1 + 



x 



n 



uniformly on any finite interval in 
\G k (X)-G(X)\ < J 



-e x ^0 
For n = 2k 

Tr(XP)\ 2/ 



(17) 



2k 



) 



_ Tr(XP) 



da(P) 



Since X ranges over a compact set, | Tr(XP)| < C for some C not depending 
on X. The integration set P is compact, so the convergence Gk — > G is 
uniform. 

Now let us prove that the sequence Bk tends to B. The minimum is strict, 
so G(X) > G(B) for any X ^ B. The sphere \\X — B\\ = e is compact, then 
there exists 

a = min G(X) - G(B) > 

||X-B||=e 
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Since G k — >■ G uniformly on compact sets, a number N e exists such that for 
any k > N e 

G(X) - | < G fc (X) < + I 

for any X : \\X — B\\ < e. Then for any such X 

\a = a-^< G{X) - G(B) - | < G fc (X) - G(S) 

So, 

+ G(B) < G k (X) WX : \\X - B\\ < e (18) 

It remains to check that \\Bk — B\\ < e for all k > N £ . Suppose a number 
k > N e exists such that the appropriate minimal point \\B^ — B\\ > e. For 
t = e/\\Bk — B\\ consider the convex combination 

E = (1 -t) ■ B + tB k 

then || £7 — B\\ = e. Then it follows from (ITS]) and the convexity of Gk that 

\ + G(B) < G k (E) < (1 - t) ■ G k (B) + tG k (BB k ) < 

< (1 - t) ■ G k (B) + tG k (B) = G k (B) < ~ + G(B) 
— contradiction. 

□ 

Theorem 8. If the sequence B k of the minimal points of G k (X) converge, 
then the function G(X) has strict minimum at point B = limi^fc. 

Proof. For any X 

G k (X) > G k (B k ) 
For any fixed X the sequence G k (X) — > G(X), so 

G(X) = \imG k (X) > BE G k (B k ) (19) 

Now check that 

lim G k (B k ) = G(B) 

The convergence Bj — > B and G k — >• G together with the continuity of each 
G k (X) imply 
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. Vft G k (Bj) ->• G k (B) 
• G k {B) 

Applying Cantor diagonal method, we choose a subsequence 5j fc such that 
Gk{Bj k ) G(B). For arbitrary X, K we have 



|G fc (X)-G fc (y)| < 



Tr(XP)\ 



2k 



2k 



Tr(YT)\ 



2 k 



Observing that for any x, y G 

X s 



1 + 



n 



- 1 + 



2A; 



< e N+|y| . u_ 



da(P) 



and taking into account that Tr(XP) < 



|G fc (X)-G fc (y)| < (( 



\\x\ 



\Y\\ 



x-y\ 
we get 

-IMD-ll-x- 



Y\ 



Using this, we have 

\G k (B k )-G k (B jh )\< (e» fl *ll+ll fl iJI 



Therefore G(B) = limG k (B k ), so it follows from (1191) that B is the minimal 
point of G(X). In accordance with the Proposition [5] this minimum is strict. 

□ 

The obtained criterion can be strengthened. 

Theorem 9. If the sequence B k of the minimal points ofG k (X) is bounded, 
then it converges, and the function G(X) has strict minimum at B = \imB k . 

Proof. Since the space £ is finite-dimensional, we can select its converging 
subsequence B kn — > B. Let us first show that B = minG(X). For any fixed 

X 

G{X) = \imG kn (X) > lim" G kn (B kn ) (20) 

Now check that 

lim G kn (B kn ) = G(B) 

Proceeding in a way similar to Theorem [6J we select a sub-subsequence B kn . 
such that Gk„(B kn ) — > G(B) and get the required: the function G(X) has 
minimum. Then, applying Theorem |5] we infer that the minimum is strict 
and using Theorem [7] we obtain that B k — > B. 

□ 
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3 "Temperature" theorem and its consequences 



The mapping 




WPda(P) 



(21) 



p 



from all self-adjoint operators to positive operators in H was shown (proposi- 
tion [5]) to be injective. However, £ is not surjective: clearly, no pure product 
state can be represented this way. We shall show that the image of £ contains 
almost all separable density operators. 

3.1 Matching theorem 

The decompositions of a given density operator p - both discrete and con- 
tinuous are known to be non-unique. The following theorem shows that any 
continuous positive decomposition of p can be replaced by an exponential 
one, which we studied before. 

Theorem 10 ("Temperature" theorem). Let p be a density matrix such that 
it can be represented in the form ([2]) 



p 

Proof. Taking into account ([7]), it suffices to prove that the function G(X) = 
J P e Tr ( xp ^da(P) — Tr(Xp) has minimum. Writing down p in the form ([2]), we 

get 



G(X) = J e Tr(xp W(P)-y w(P) Tr(XP)da(P) = J (e Tr(xp) - w(P) Tr(XP)) da(P) 
p p p 

(22) 
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For every fixed P 6 P introduce a function 

g {t) = e*-at (23) 

where 

a = w(P) , t = Tr(XP) 

The set P is compact, and the function w(P) is continuous and positive, that 
is why it attains its extreme values 

< m < w(P) < M 

Since w(P) is a probability density, we have 

< m < 1 < M (24) 

The value of a in (|2"3"|) lies between m and M. The minimal value of g(t) 
depends on the value of the parameter a as follows: 

fiWn = a — a In a 

By elementary (for t < 0) and routine (for t > 0) calculations we obtain 

g(t) > (m + M) [1 - ln(m + M )] + m|t| (25) 
for any real t as illustrated at this graph: 



50- 
40- 
30- 
20- 

-3 -2 ' ' 1 ' 3 






~* 555 ^^ZIx^Iai^ 



Since t = Tr(XP) and a = w(P), this means 

e^ 15 ) - w(P) Tr(XP) > (m + M)[l - ln(m + M)\ + m\ Tr(XP)| 
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Integrating this inequality over P, we obtain the following evaluation: 

G(X) >(m + M)[l — ln(m + M)]+m J | Tr(XP) |d(x(P) (26) 

Consider the second summand of this expression: 

u(X) = J |Tr(XP)|da(P) 

which is a seminorm. Proceeding like in Proposition [6] we see that u(X) is 
non-degenerate. Since all norms in finite-dimensional space are equivalent, 
u(X) evaluates as 

v[X) > const • \\X\\ 

So far, we obtain 

G(X) > (m + M)[l - ln(m + M)] + const ■ \\X\\ +00 

Since G(X) is convex, it has minimum, which is attained at certain point 
B e£. □ 



Remark. It can be verified that the distribution w (P) = e Tr( - BP - ) yields the 
maximum for 




w In wda(P) 



Theorem 11 (Matching theorem). Any p belonging to the interior T> of the 
set of separable states can be represented in the form (j3J); 

p = J e Tr ^ p )pda(P) 
p 

Proof. Denote by Ai the set of density operators representable as Husimi 
exponentials. Let pi, p 2 G Ai and form their convex combination p — (1 — 
t)pi + tp 2 . This p can be represented by continuous positive density w(P) = 
(1— t)e Tr ( BlF ^ +te It ( B ' 2P \ therefore it follows from the "Temperature" theorem 
[TD] that such B exists that p = J e Tr ( BP )Pda(P). That means, M. is convex 
together with its closure. 
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Let Po G P be a product pure state. Any atomic measure on a compact 
set P can be approximated by a sequence of strictly positive densities. Take 
such a sequence w n on P such that for any continuous function / on P 



lim / /(PK(P) da(P) = /(P ) 



71-> 

P 



Consider the sequence of operators 

w n (P)Pda(P) 



each belonging to M.. For any affine function h(X) 

h(p n ) = J w n (P) h(P)da(P) -> h(P ) 



p 



therefore p n — > Pq. So far, the closure of Ai contains P. The closure is 
shown to be convex, therefore it contains all separable density operators. 

Now consider the "partition function" 

Z(X) = J e Tr(xp W(P) (27) 



Note that for any 1/^0 

d 2 Z(X + Vt) 



WV) - 2 



e' n '< xp >(Tr(l/P)) 2 d CT (P)>0 



t=0 

p 



whose positivity was established in (1121) . However d 2 Z = dC, where C = 
VZ, see (T2T]). Then the mapping C(X) is non-degenerate at every point 
X G S. According to the inverse image theorem, we see that C is local 
diffeomorphism, that is, each point X has a neighborhood whose image under 
the mapping £ is open in E. So, we conclude that its full image £(£) is an 
open subset in £. In the meantime 

M = C(S) D {X G £ | TrX = 1} 

so, M. is an open subset of the set of separable density operators: M CD. 

So far, we have shown that the set M. is an open, convex and dense subset 
of the set of all separable density operators, therefore M. coincides with its 
interior: Ai = T>. □ 
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3.2 Husimi coordinatization 

Starting from the "partition function" Z defined above (1271) . introduce 



W = \nZ = In J e Tr(xp W(P) (28) 
p 

and first calculate its gradient 

W (X) = ^ = I J e^^Pd^P) 

p 

Proposition 12. T7ie quadratic form d 2 W is nonnegatively defined and van- 
ishes only on scalar operators. 

Proof. Calculate the value of d 2 W at point X on element V 

J_ />«*+«>P>Tr(l/P)d<7(P) 



Denoting 
it reads 



- I 

dt Z 



t=o 



Tr(XP) 

w(P) = —— (29) 



d 2 W(X;V) = J (Tr(VP)) 2 w(P)d ( r(P) - (j Tr(VP)w(P)d ( r(P) j > 

(30) 

This expression vanishes if and only if Tr(VP) = const for every P 6 P. All 
such V are the scalar operators: V — XI. □ 

Note that for any scalar operator the function W has the property: 

W(X + XI) = W(X) + X 

therefore its gradient 

X (X) = VW{X) (31) 
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is invariant with respect to shifts along scalar operators: 

X(X + \I) = X (X) (32) 

This shows that the mapping x is n °t injective on £, but its restriction to 
the set M of traceless operators from £ becomes injective. 

Theorem 13. The mapping fl3~Tj) establishes a diffeomorphism 

X:M^V (33) 

between the set M of all traceless self-adjoint operators in H and the interior 
D of the set of all separable density operators. 

Proof. The mapping x was shown to be invariant with respect to shifts along 
scalar operators, so x(AT) = x(£). The differential dx = d 2 W . It follows 
from Proposition [12] that d 2 W > on A/", therefore dx is non-degenerate 
on M . So, x is a diffeomorphism of Af on the image x{N). In turn, the 
image of x coincides with the image of the mapping (1211) - - the density 
operators representable as (jl]). Finally, as shown in Matching theorem [TTj 
this set coincides with T> — the interior of the set of all separable density 
operators. □ 



Concluding remarks 

In this paper, the standard thermodynamical ideas were employed. However, 
it was not possible to replant directly the standard techniques. The first 
reason is that the analog of inverse temperature is no longer a scalar 0, it 
becomes self-adjoint operator B. That is why the "Temperature theorem" 
needed to be reproved in the new setting. One more essential difference is the 
behavior of the "partition function" Z. Unlike classical case, in our setting 
the gradient of \nZ is invariant (13 2 p with respect to shifts: X \- > X + AI. 

The overall construction is general enough: we do not dwell on any par- 
ticular decomposition of the state space H into tensor product — this notion 
was shown to be relative to measurements [2]. Moreover, the set of pure 
product states P could be any compact full set of unit vectors in the state 
space H. 'Full' means that the value of any quadratic form in H is com- 
pletely defined by its values on P. This is the central point of the method, 
based on multipolarization identity f fTUj) . which was inspired by the ideas of 
Kodi Husimi [3j to consider a special kind of functions defined on pure states, 
namely, those parametrized by values of appropriate quadratic form. 
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